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LT eODUOTION 


The role of the subcritical assembly is gradually being 
enlarged to include studies of transient or dynamic resnonses. 
Demands for an inexpensive and quick rethod for determining 
the offect of various paranetears upon control reszonse have 
led investigations into the utilization ef the assemoly Lor 
guch purposes. Perheansa the most important advantaye of the 
suberitical assembly is the ease with which exverinsnts may 
be performed. Uxtra physical safeguards and restrictions are 
at a minimum, and no extensive shielding is resuired. The 
relatively small size of the suberitical assembly also tends 
to make it easier and leas exvrensive to vary the desired 
parame ters. 

Borrowing from servomechaniam theory and treating the 
assembly as a “black box", many investigators have aptly 
applied the concent of the transfer function to describe the 
response to a given forcing function. In practice, phase and 
amplitude relations, whon considered as functions of the 
driving frecuency, are collectively designated as the assemoly 
transfer Punction. 

The parametric forcing or driving function is of two 


types, reactivity foreing and source forcing. In general, 
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feaativity foroing 18& eepresented by the chanse in neutron 
gensity due to a small oarmonic variation in the reactivity, 
that ic, go harmonic change in the renetor geometry accon- 
piished vy oseilluting some neutron absorber. Jouree 
forcing, in contrast to resetivity forcing, in no way affects 
or chenges the assembly geometry. it 4s merely an imposed 
V@iatiou in the primary sourve. 

The objective of this thesis is to present the sub- 
eritical anserbly transfer function which was determined by 


tie technique of source forcing. 





iMVISw® OF THe LITMUTUAS 


The history of the subcritical pile and its role in the 
Gevelopment of the first self-sustaining chain reaction has 
become legend (1), but work is continuing (2) to make the 
agsendly a testing vehicle and not merely an object of in- 
vastication for itself. 

Tho initial servomechanism concest was presented by 
Figott, Crever, and “wens (9 and 10) at the Knolls Atomic 
Power Laboratory. Because the roots ef the kinetica enuation 
are negative, the vile tranafer function was noted ag veiling 
the sane as that for a minimum phaee shift network. This 
permitted it to be rapresented by a Bode diagram. aA pile 
kinetic simulator, shane voltace rerponze represented the 
pile neutron fiux responce, was successfully utilised for 
me test to determine the stability of an automatic control 
System and ites ability to maintein a desired power level. 

Franz (5) determined susstantially the sane tyne trans- 
fer function but used a more sephisticated approach. By 
definition, the transfer fumetion relates sinusoidal varia- 
tiens; hence it can oe handled by the standard techniques of 
servomechanion theory. The ratio of the Laplace transform of 
the output to that of the input waa found to give the transfer 


funetion more conveniently. in both the above situations it 
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was oarefully noted that the functiona determined were those 
ef a pile simulator and described the pile reaponses only to 
the same degree that the ecuations utilised descrived the 
conditions in the pile. 

It remained for Harrer, Boyar, end Rrucoff (8), based on 
their now classical experiment on the Argonne C?r-2 reactor, 
to justify the use of the transfer function to describe a re- 
ector as mn control leon element. ‘Since their experiment it 
has become standard procedure to determine the transfer func- 
tion of critical assemblies and reactors for stasllity eri- 
teria. This method has also been applied to mockups of pro- 
posed reactors. Soland, Gnith, and Kice (3) found the zero 
power transfer function of the uPi-<III mockup assemoly of tire 
WikeI, a fast critical assemly. The results were used to 
determine the best oscillator roda and nreasurement techninues 
to be used on the actual assembly, when constructed. 

In all the veferences found, the information primarily 
pertained to a critical assemoly and to the solution of the 
transfer function by reactivity forcing. OUtner than a cur- 
sory discussion in Glasstone (6, p. 252), no unclassified 
reference was found which contained any information con- 
eerning the transfer function of a sudcritical aasembly as 


such or concerning souree forcing. 
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Dee wVal Toe OF WEE TRANS Pi) PURCTION 


The response of the suboeritiscal unsembly to scurce 
forcing cam be derived from the time-dependent general diffuse 
sion ejuation. For simplicity all neutrons are considered 
thermal, end hence the e mwiation can be written as (7, p. 


201) 


db g@t- 2 g+s= Sh - $e Lge 1 


where nr is the meen velonlity of the thermal neutrons. in 
general the symbole in the derivation are in common usage (6, 
PDs 223-253), but the moanings are given in the List of 
Symbols for ceomrrentont reference. 

The source term is made up of theese rarts. The Sirat 
portion ia the sromet neutrons resulting Prom Pisaion. 
Aeoording toa she Fermi are theory thia promt neutren source 
tome nay be aexpregsed ag 


> “n't 
S,*° =f) k, a %e Ba. 2 


where (3 2 the total fraction of delayed nmitrons in any 
one generation. The seceord portion is the delayed neutrons 
which result from flasion and which are given of f during the 
decay of certain fission precducts. Generally, six groupe are 


recognized, @¢ach growo having a dinerote ceeny constant r, 
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whist of dynbois 


wuckling om, 7¢ 
Praotion of average group of delayed neutrons 
Yraction of 46h group of delayed neutrons 


Concentration of avernge group of delayed 
neutron precursors 


Goneentration of 1>h group of delayed neutron 
PYeCursgel s 


Thermal neutron diffusion coefficient OMe 
bffactive multiplication factor 

Multiplication factor for an infinite assembly 
Reutron diffusion iength in the moderator CMe 
Average neutron Lifetime in a finite assembly sec. 
Mean neutron Lifetime in on infinite assembly sec. 


Decay constant of averages group of delayed 
neutren precursors sec .7t 


Decay vgonetent of the 4h group of delayed 
neutron precursors sec. 7+ 


Mumoer of neutrons of thermal energies 


hesonancs escape propability 


Thermal neutron flux NOU trons 
om.” see. 
kate of production ef source neutrons neutrons 
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Macroscopia absorpti.n crosssectian emt 
Time SEC. 
Fermi ace or. © 
Thermal neutron mean velocity em. /SeC. 
Fremuency of oscillation radians 
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and fraction (3, (11, p. Uy). With C, defined as the concen- 
tration of the 1th group of delayed neutron emitters, the net 


rate of formation of the precurnor in given by 
Kia 


Since the rate of eS of the delayed neutrons is 
numerically ony at to She decay rate of the precursor, thie 
rate ecuals 2 Ay Cy for the six groups under considera- 
tion. itn the ermi «age of a dolayed neutron assumed equal 
to that of a promot neutron, the slowing-down nonleakage 
probability will be _ Bet This is net strictly true tecsause 
the energy of the delayed neutrons is senerally less than the 
average energy of the prompt nevwtrons, but the approximation 


is considered reasonable (6, p. 227). The delayed neutron 


source term may be written as 


2 iS 
ct “pT oa 


whe re 7 is defined as the resenunte escape probability. ‘The 
third portion is tre external source term which will be 
defined as 
s+ 3,+8,0°°* Eq. 5 
Wnere 5, ‘is the complex amplitude. 
Yor further simplification of the mathematica, the apace 
aid time variables will se consicered seperable so that the 


partial differential equations may be reduced to ordinary 
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differential equations (6, p. 227). 
With substitution of Bquations 2, 4, and 5 into wcquation 


1, the resulting expression is 


6 
2 2 
2 = CU - me 
DY - aa + (le Qik, 2 je . + Pos ; Pa Ao, 
4 (S.+5, gi@ ty - Fr se. qe & 


Only first modes will be considered for the simplified solu- 
tion, and hence 7 >= -3% (8, p. 33). By definition ‘Le is 
the mean lifetime of tue tnermal neutrons in an infinite 

es - With division by 2. and 


a 
the above substitutions, ouation 6 becomes 


medium, and it ecuals 
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where k is the effective multinlication factor and y is the 
averaze thermal neutron lifetime, both in a finite assenoly. 
With term rearrangement and with substitution of k and AL P 


1t is possible to weilte Equation 7 as 
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d 
with division vy L » collection of the terms containing bs 


and substitution of @ = nar, vhere ar is considered constant, 


Equation 10 reduces to 
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With these substitutions and division by/fu", HKquation 11 


yecomes 
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With rearrangement of Bauation 15 and with substitution for 


Agfa, i quation lt besomes 
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with source forcing, k will be a constant ocecause the 
bile geometry is not chanred. The corressonding forcing 
value of neutron density will be 
> ie oye’ * BQ. 19 
were Ny is the steady component and Ry is the complex ampli- 
tude. in an analogous manner 
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By combination of Bguations 15, 19, and 20, Mguation 18 
veoomes 


6 
2 (3... 
NM, + nye = 1 i ae (n +n,@ ) 


(k= 1) ~ 
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Only one average grouo of delayed neutrons will be con- 
sidered in order to maze the equation more wieldy. This 
introduces some additional, dsut small, error inte the results 
(6, p. 251). Houation 21 is the complete form, but by the 
definition of a transfer function only the sinuscidal portion 
is of consern. With separation of the steady and sinusoidal 
ecemponents and with the srou>p notation dropped, the 


sinusoidal vortion of Equation 21 5ecomes 
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The simisoidal portion of Lquation 23 becones 
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With the transfer function defined as the ratio of the 
response to the sinusoidally varying driving foree, the 


desired ratio is 
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DRVORIPTIUON OF ah Pa ATUS aX) BAPHAIMLNTAL ZeMOCEDUUS 


Figure 1 is an overall view of the subcritical assembly 
with the east cover removed. iteactor grade graphite 
oylinders 60 inches long and 7 inches in diameter were 
machined to a square form 6 inches across, leaving rounded 
corners. These were stacked in ten columns nine rows high. 
Additional cylinders were eut to a S-ineh high by 6-inch wide 
form for the top five rews. The final rile sige, with the 
machined cylinders stacxed fourteen rows high and ten colunns 
wide, measures 60 by 60 by 79 inenes. The whole array of 
graphite is surrounded by a sandwich of plywood, 0.010 inch 
eadmium, and masonite, bolited in sections and braced with a 
fremework. The east and west covers are removable to gain 
access to the pile. 

The graphite is placed on a wooden base constructed so 
that there are three accersible spaces under the pile. These 
Spaces extend the north-south width of the assezbly. Two 
aluminum tanks, also extending the width ef the pile, were 
filled with water and placed in the two outermost spaces to 
moderate the neutrons and to reduce radiations coming di- 
reetly from the sources beneath the pile. The central apace 


was used to house the oscilletor shown in Figure 2. 
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Figure 1. The sudcritical assenbly 
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Figure 2. ‘The oseillator assembly 
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The holes between the graphite blocks, measuring 
approximately le inches across, permit insertion of the 
fuel elements in varicus lattice arrangements. Aluminum 
cylinders approximately 62 Inches long each containing seven 
8.2-ineh uranium slugs constitute the fuel elements. 

The east cover was in place during the test runs, and 
five holes were drilled on the cover centerline at distances 
of one foot, two feet, three feet, four feet, and Live Leet 
frem the base. For all further discussion these holes will 
be designated by ordinals with the hole one foot above the 
vase veing called the first hole, the hole two feet above 
the base being called the second hole, and go on. With the 
8.5-inch lattice, as shown in Figure 1, the cover holes were 
in line with eomty holes between the graphite blocks, per- 
mitting a neutron detector to be inserted into the pile, 

The apparatus used was assembled from readily available 
material and parts. xtreme aceuracy and refinenent of 
Gesign were not attempted. Five Puebe sources cylindrical 
in shape, each measuring 1 inch in diameter and 1 3/8 inches 
high, were placed in the source retainer of the oscillator 
shown in Figures 2 and 3. Moderation was provided by a 
paraffin jacket formed about the source retainer. The 
jacket was made by rolling 1/16-inoh aluminum sheet to an 
outside diameter of 7 3/8 inches and cementing it to 1/3-ineh 


end plates. This container was filled by means of a filler 
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hole with 25 peunds of high melting point white paraffin and 
slowly cooled to minimize the possibility of shrinkage voids. 
Approximately three inches of paraffin swrerounced the sources. 
The retainer acted as support for both the non-rotating 
paraffin jacket and the bearings on which the eaternal 
eylinder rotated. 

The rotating cylinder wes fermod by rolling 1/8-inch 
aluminum sheet to an cutside dlareter of @ inches and cutting 
it to a length of 15 inches. The end plates, also of 
aluminum, were from ,; inch stock, shrink fitted and screwed 
to a collar fitted on the bearing outer run casing. ‘The 
cylinder was screwed onto the end plates permitting dis- 
assembly. The sine pattern was cut from a 0.010-inch thick 
sheet of cadmium to the dimensions shown in Pigure 3 and ce- 
mented with an aluminum vased cement to the inside of the 
rotating cylinder. 

Por further discussion position 1 will be used to refer 
to the rotating cylinder and pattern pesition whenever point 
A of tne flattened pattern in fipure 3 is top dead center. 
Pesition 2 will be used to indicate that point bl is top dead 
center. Unsed on the difference of readings taken with the 
cadmium pattern in the two extreme positions, it is estimated 
thet approximately 20 per cent of the source neutrons reached 
thermal speeds in the paraffin. It is also eatimated from 


the accepted macroscopic absorstion cress section of cadmium 


“~ 


due ACNSC eed othte omla At” ad fo cues Se eer seme 
thle QA EF) YPN ee ee mete ot fer eee 
eT AO mae ALT Ce met Yim 
Oa 
ee ee ee ee eT | 

i ome * ne . me bry tls ge 
getline Cie eo eee ee 
paler oth fe dt ah ot Ret ea 
ob eee ee or ee YS ome ort 
pees vie edi) dee ee et ey? vee opel 
+ winter were 5 mur new OW 














20 


ROTATING CYLINDER SS 


WITH PATTERN a 


SOURCE RETAINER 












NONROTATING 
PARAFFIN 
JACKET 


WOOD EXTENSION 


POINT B 24. 35 IN. 


Figure 3. Oscillator and cadmium pattern detail 
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that 93 per cent of the thermal neutrons reaching the pattern 
were absorbed, 

kn extension of one of the bearing collars was ;rooved 
to hold the leather driving belt. The rotation frequency was 
varied by eontrolling the voltace to a +-inch Black and Decker 
electrics drill mounted on a four-foot wooden extension. A 
Speed reducer was interposed to permit speeds down to approxi- 
mately 5 RPM. ‘The oscillator was desirzned for speeds an high 
as 1060 RPM, but no tests were carried out in this rance. 

For maximum signal the neutron detector, a 3/3-insh 
dianeter, Seinch long tube filled with enriched BF, gas, was 
inserted in the first hole to the pile midpoint. A Radiation 
Instrument Vavelooment Laboratories, model 200eI, scaler was 
used for the high voltaze and initial amplification. The 
amplifier-discriminator jumper was disconnected and the 
amplifier output taken directly to a Nuclear Chicago model 
1615 survey meter with the high voltage rendered inoperative. 
The output of the survey meter was fed to a type 112 direct 
couplec Tektronix amplifier and thence to one amplifier and 
pen of a two pen Srush recording system. The instrument 
arrangement is pictured in Figure lL, and all the circuitry 
utLlized is shown by block dissram in Figure 5. 

The time applied variation in source strength was re- 
corded by means of brush contacts which were mounted on the 


oscillator frame and the rotor and which were wired parallel 
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Figure 5. Electronic circuitry 
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to a 1lOK resistor as shown in “igure 5. The 10K resiator 
was in series with a 1004 resistor, and a 110-velt 40- 
cyele voltage was impressed across them. The second re- 
corder amplifier and pen were also connected across the 10k 
resistor. With the rotor contacts open, a 1l0-volt 60-cycle 
fmput signal was plotted by the pen. whenever the rotor 
contacts closed, this signal was shorted out. The orush 
contact mounted on the frame was approximately 120 dcgrees 
from rotor top dead center wnile the rotating contact was 
set shout st point A on the rotor. It was found that 
oscillator speeds up to about 300 KPH could be measured in 
thig manner. At faster speeds the length of time tho wave 
{es shorted is only a portion of one wave cycle and is very 
difficult to read with the rather slow pen response, 

neadings were taken with the oscillator steady in posi- 
tions 1 and 2 to determine the maximum range. This simulated 
an infinitely slow oscillation. The difference in the 
readings was called 5,'. .eadings were then taken with the 
osclilator rotating at various speeds, and for each fre~ 
quency the peak to peak amplitude, called n,', was compared 
with the base value 5,'. 

The phase angle was determined by plotting position 1 
of the rotor on the time wave and comparing it with the 
position of peak flux for each cyole. It was noticed during 


preliminary tests that moving the orobesteadily in and out 
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of the pile by hand produced a recorded wave having an 
amplitude les# tnan that ultimately recorded if the prope 
were permitted te stop in the pile for elifnt to sen aeconds. 
The pen, however, seemed to resvord in pnase with the motion. 
It was therefcre decited to correct the recorded amplitudes 
for the two second time constant of the survey meber and to 
accept the recorded phase angle as veing aporoximately 
correct. 

feeanuse tne theoretical cireuit analysis appeared to te 
too lengthy, it was decided to ovtein an amplitude correction 
experimentally. a sine wave signal generator was connected 
to the survey meter input and adjusted to give a e2=-volt out- 
put. The scaler amplifier was not included in the test cir- 
euit beceuse of the internal elreuitry modifications that 
would have been entailed with t..e inclusion. since tne 
acaler arplifier was in essence amplifying spike voltages 
with a five miecroesacend resolution, it wae considered to have 
negligible effect on the variation of the recorded signal 
amplitude. ‘The signal generator frequency was inereased 
until the survey meter registered 20,000 counts per minute 
whieh was approximately the same as the reading obtained with 
the probe inserted to the pile midpoint in the first hole. 
It had been noted prior that as the oscillator was rotated 
very slowly, the survey meter reading varied between approxi- 


mately 19,000 counts per minute and 21,006 eounts per minute. 
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Tho aignai generator was adjusted until the survey meter 

read these vaiues, and the extremes were marced on the dial. 
The sinusoidal variation voetween these extremes was ovltalined 
by atteching a yardstick solidly to the vignal generator 

dial to vermit a pendulum action. aA heavy weight wes cisamped 
at various distances on the yardstick and swung to simulate 
the different cycling frequencies. With the pendulum held in 
the extreme positions the amplitude difference, as sefore, 
Simulates an infinitely slow cycling sveed. The ratio of the 
observed peak to psak amplitude at any one frequency to the 
aaplitude difference ceterained above was talcen as the oor- 
reotion facto: for that fremency. The arplitude correction 
curve is plotted in Figure 6. The rmee of readingn for the 
correction eurve vas sevorely limited by the pendulum length. 
Because Of the physical arrangement, pendulum lenetas lers 
then 15 inches were not feanible. This was not a serious 
handicap sinee shorter vendulum langths cerraspended to 
faster eycling soseds which were beyond the upper values of 
the experimental data. Ths 30 inch limitation of the yard- 
stick, caured by mounting on the dlal, corrgssponded te a 
frequency of 3.59 radians per second using simple pendulum 
theery. It is to be noted that to obtain a freguency of one 


recian per second « 32.2-foot pendulur would have been 


needed. The correction wae applied by Gividing n,!' by the 
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AMPLITUDE CORRECTION FACTOR 


Figure 6. 
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Amplitude ratio correction factor vs. frequency 
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corresponding correction fasteo: at the same frequency to give 
the experimental values of amplitude patlo, n,/S,'s plotted 
in Pigure 7. 

Some test mins were attempted using the sesond and third 
holes but the data were not readable. Apparently the 
Statintical variations were of the sano order of magnitude 


as the flux variations, masking the desired signal. 
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DISCUSSION orf  2ASULTS 


The primary results of the experiment are presented 
grephically in Figures 7 and %. The data from which these 
figures were drawn are tabulated in the Appendix. 

The smooth curves are plots of Equation 25 for different 
values of the mean neutron lifetime in a finite assembly and 
for different values of k. Since neitiner L one k have been 
determined for the assembly in question, values of k = 0.7 
and 0.8 and X = 0.07 secormis and 0.1 seconds were chosen to 
indicate the trend of the curvas with both yi ana k ard to 
indicate the general curve shape. In the theoretical deriva- 
tion the six groups of delayed neutrons were replaced with an 
everege group to reduce the amount of computation. With the 


total faction (3 of the delayed neutrons equal to 0.0073, an 


averace dr equal to 0,08 seconds” * was computed from 
2; 
= % Piri (6, p. 230). With the above paraneters 
3 


in the Equation 26, the amplitude ratio and phase angle were 
determined for verious frequencies. 

Measurements of the neutron flux in the first hole were 
taken for oscillator speeds from approximately 0.5 to 5 
radians per second. Both ohase angle and amplitude ratio 


readings were obtained. The secorded phase angle is plotted 
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Figure 7. Amplitude ratio vs. frequency 





jin Figure 8. The recorded amplitude ratios were corrected 
for the two second time constant of the survey meter by the 
correction faecter indicated in Figure 6. The data for the 
amplitude ratio sorrection factor, tabulated in the Apvendix, 
plots as a straight line on semi-logarithmic paper. The 
range was increased by straight line extrapolation. This 
result is valid oniy over a limited range since for all fre- 
quencies the factor cannot be greater than unity or less than 
zere. lence the true factor must deviate somewnat from a 
atraignt line as it approaches the extreme values. at the 
lower frequeneies the correction factor itself becomes small, 
so the abeve deviation will have little effect on the cor- 
reoted amplitude ratic. At the higher frequencies where the 
correction factor becores large in cosiparison to the readings 
the effect in magnified, and extreme care would need to be 
taken to justify any extrapolated values. Fortunately tie 
frecucreies obtained in the experimental correction factor 
tests wore as high as those used in the amplitude ratio tests 
and no extrapolated values in the high frequency range had to 
be used. The w plotted in Figure 6 is actually that of the 
pendulum, but it corresponds directly to that of the 
oscillator since in both cases it ie the frequency of the 
Griving function. The fremency was solved for a simple 
pendulum (without damning) from the formula Ww = fey 


were g is the gravitationnl acceleration in feet per secon d® 
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Figure 8. Phase angle vs. frequency 
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and p 48 the pendulum length in fee re 

AN interesting result of the experinent was the decrease 
in median flux level with increasing GW. The value of the 
median line of the sinusoidal flux variation for each testing 
frequency was plotted in Figure 9 as it appeared on the pen 
recorder, with the recorder paper scale ag the arbitrary 
ordinate scale. ‘This value corresponds to Ry in the theo- 
rotical development, and it is tabulated in the Appendix for 
each frequency tested. In the splitting of Equation 21 into 
the steady and the sinusoidal conditions, it is seen that De 
has no theoretical relation with ©, lence tne odoservation 
wags not predicted by the theory. The only legical explana- 
tion lies in the oseillator design and the method chosen to 
accomplish source forcing. No tests were run to dotermine 
the exact ecmize since the actual median flux level would have 
little effect on the cycling amplitude ratio. 

Experimental uncertainty results from many factors. A 
very slow rotation trial run was made with the neutron probe 
within three inches of the oscillator to check the variation 
3n the source. The very small distyurtions noted in tne 
plotted sine wave could possidly have veen due to neutron 
scattering, voids in the paraffin, and amplifier drift. 
Perhaps ths greatest source of error in the experiment was 
tne signal fluctuations due to statistical variations and 


noise. These variations obscured an) distorted the peaka of 


aed i eet ett 6 ot & On 

OO ee ee 
=i hm whey oe - 6S entacemund pisy teond wal mulvee af 
en i a 
Ge ee ee i ee 
Teel Me od bb chee ee eee 6 Oe veer 
rae onl 6b ff Ot tlmmives Milas Gime ee es ae 
ee ee ee ee ee 
AAS IS eelieeet Te pale ie) oe eh OEP NY arene ry Ge 
ff Rett Ae AP AL panned eco Legis me Hew qeenne oad 
everest ee ee a ee 
saute deete! gee ee? ogee eee oe ee 
9 tes ume ot Lm cg ico tame i liane a@ a! erdh ale | | 
OS At me Rete ee rte hd hom - 
ce blew deems cal lee Lee ow ome omer Seve 
«other! tvibens sues ae ay BelTS oc208t 

— comme gee wwe) me teeny pI ede | at ee 
ee ee ee ee 
aebdetGee ott fee Of SEttia a S ete @ee ola 
ee ee ee ee 
ents 70 1 cae! Shed Pee Ghee Gee eee Ceeee 
EVIE TATRLL Om. ow Yeh Ttrmn ee wd calee guelverdone 
RAP Reine we) OT Rat he Dee et om 
he MOLAR NAAR LEG RR OP Gm mtn arenes Comzee 
“2 Get ee ate Pate a tee been Genta ae met 
































ARBITRARY SCALE 


MD 


nn 


f 


OJ 


04 


O6 08 l 


35 


Z 4 6 8 lO 


FREQUENCY - RADIANS / SECOND 


Figure 9. 


Median flux level vs. frequency 


20 








ee er eel eel 





Pigure 10, Yiews of typical data 
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the signals, adding to the difficulty of determining both the 
gwaplitude patio and the phase angle. Two samples of typisal 
dute, demonstrating readings obtained with both high and low 
gain, are shown in Figure 10, 

To compare the experimental and calculated values of 
aeplitude ratic and phase angie, it is only necessary to 
observe the sraves of the curves. it is seen strat the 
exporimantal peinta lie reasonably elose to the theoretical 


curve plotited for values of « = 0.8 and XL = 90.07 seconds. 


= wee etnies > eT a toe ee ae od 
Maton 1 rey ote GS me Slee vil 
wl Gs Gat ow ND Beer atti . Ante ote 
2 ee! a) oe we ey 

I ee el 
eat cee I ele ee oe eel 
Oe ee ee 
0 od eel hs ee a de A 
dew Gel + ¢ We Games Hy baAaie 




























Yo 
OD 


ORC LUSIONS 


The transfer function determined experimentally agrees 
quite well with that predicted by the theory. hile it was 
not the purpose of the thesis, the experimental points indi- 
cate a mean neutron lifetime in the assembly of approximately 
0.07 seconds and an effoctive multiplication factor of 
approximately 0.8, 

The roter assemoly thermalized and absorbed a sufficient 
percentace of neutrons to permit obtaining acceptable data 
from the first hole. In the second and higher holes the 
Sinusoidal variation in the neutron flux was reduced to the 
order of magnitude of the statistical and noise variations 
and could not be read. With a larger or more efficient 
neutron deteetor, a thicker paraffin jacket, or possibly a 
thicker cadmium pattern, the amplitude of tha sinusoidal 
variations could 5e substantially increased permitting 
correspondingly smaller deviations and statistical errors. 

The median level of the pile neutron flux decreased 
with increasing cycling speed. This result is not predicted 
by the theory and ig p.obsably die to the desizn of the 
escillator,. | 


Althougn the experimental results were not extremely 


— - 
Sore LAemmmins Meatindeh hums nthe aA 
om hte eptte ee 14 edidhoy, Yue tof? Line ety 
thind eanlens Aalomnl vets 28? yritanr 8) o> Keyan owt oe 
Mel aahety Te \etmane 4f% 41 AMAAET) som teee eee © eee 
a RE aa ALO wT we ome come UA 
hig (hte 
ORR T et 6 frre cee beckett 5 eee veda aol? 
meet elimtyeme tndties Fires ee fee eb 
<a) eee woul os ces oe) met fm a aS 
rd bg bene ee eT) ee ae tet me tom 
saedeminse eohts me beeh/dicnde Ot Te abinimar Se aame 
dmetAiTe oon @ ee 6 Oe oe ee eee bee 
0 ieee © dee A/S I ert eee eee 
Ioheweende oe SE we ehom ool gm Anns omrmen comates 
tO me meres i LLAmm ate tee we heme wmntenseee 
eewrers (Me Mal Ghdy aoe CENTS) FY CORR fine ree 
remover: oolT crue wlig we Te ted mite of 
sk re 1 freee sie com dine: atte unl ate 
a & Glee bet 8s eee (Oke © OF ee pet ae Ww 


ee 





ee 





oO 
“tee 4b 1s Colt that tho data ean be reproduced with 
fPictont aecuraey and consistency to substantiate the 






one ~ 
eo ee 





Wi 


SUGORST IONS PU Pol TIE ATUDY 


It is very apparent that the possibilitiss in cycling 
@ subcritical assembly have not been fully expiolted. It re- 
mains to improve on and enlarge the range of readings to much 
faster and much slower cycling speeds. This will roquire 
some modification in the oscillator drive mechaniam. 

Neadings with a higher degree of certainty at the greater 
speeds can bo cotained oy redesigning the vasic circuitry to 
have a lower time constant and perhaps to include a dis- 
eriminator to veduce the statictical varistion. vUnce this 
nas been accomplishcd it woule be desirable to obtain tne 
transfer function of the sare pile by the techniqe of re- 
activity fercing and to compare the results of the two 
methods. 

All tne suggestions above certain to the assembly 
iteaelf as being the object of investigation. With the 
backround knowledge of the transfer function, the whole 
aroa of investigating the utiligation of the subcritical 
assenbiy as a facility for testinz components and fer 


sdvanced design researen will heave been opened. 
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Yable 1. amplitude ratio correction factors 








1 1 as Gd) Correction  -orrection 
inches feet (= ¢/1) redians/second amplitude factor 





30 2.5 12.57 3.59 0.95 0.19 
25 2.08 15.5 3.94 0.75 6.15 
20 1666 19435 el 0.6 0.12 
15 16.25 2557 5,08 0.35 0.07 





Teble 2. Amplitude ratio data 








- oe. a ah a 
1 1 
12 Oah9 0.85 0.85 0.95 
13 0,87 0.65 0.65 0.92 
Ly 1.26 0.5 0.5 0,88 
15 1.73 O.! 0.4 0.89 
16 2025 0.3 0.3 0,82 
17 2.78 0.2 0.2 0.69 
18 3.31 0.15 0.15 0.65 


19 3.73 O21 0.1 0.59 
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Table 3. Phase angle date 





Run Ww Phase angle-radiens 
radians/second 
1 Oo lpy 0.07 
2 0.715 0.27 
3 1.15 O oily 
b 1.73 0.53 
5 1.85 0.62 
6 3.31 0.88 





Yable 4. Median flux levels 


GILL: AO LL ED A Ly PS IP DLE RS EE GS i a I i a EE EE EE ee ee ke 


WW Level 
radians/ second 

Weld 73 

0,87 To3_ 
1.26 Tee 
>i Tea 
2e7 Tea 

3631 7el5 

3.73 7e3 

; z 

9.67 7.0 

12.2 720 
13.9 629 
16.h 6.6 
17.9 6.8 
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